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Abstract 

We investigate the doping dependent low energy, low temperature (T = 5 K) properties of nodal 
quasiparticles in the d-wave superconductor Bi2.iSri.9CaCu208+<5 (Bi2212). By utilizing ultra- 
high resolution laser-excited angle-resolved photoemission spectroscopy, we obtain precise band 
dispersions near Ep, mean free paths and scattering rates (F) of quasiparticles. For optimally 
and overdoped, we obtain very sharp quasiparticle peaks of 8 meV and 6 meV full-width at half- 
maximum, respectively, in accord with terahertz conductivity. For all doping levels, we find the 
energy-dependence of F ~ \lo\, while F(a; = 0) shows a monotonic increase from overdoping to 
underdoping. The doping dependence suggests the role of electronic inhomogeneity on the nodal 
quasiparticle scattering at low temperature (5 K < O.OTTc), pronounced in the underdoped region. 

PACS numbers: 74.72.Hs, 74.25.Jb, 79.60.-1 



1 



Extensive studies over the past two decades have indicated a d^^-yi symmetry of the 
superconducting gap in the Bi-based high-Tc cuprates.^'^ The dj;2_j^2-wave superconducting 
state is characterized by hne nodes along the (±7r/2, ±7r/2) direction, and allow the exis- 
tence of quasiparticles crossing the Fermi level even in the gapped superconducting ground 
state. Early in the literature, theoretical studies have discussed the behavior of nodal quasi- 
particles in terms of massless Fermions with a Dirac-cone dispersion.^ By investigating such 
nodal quasiparticles which are responsible for the low energy excitations, one can acquire 
insight into the characteristic features of the scattering and interaction mechanisms remain- 
ing in the superconducting state. Experimentally, the properties of nodal quasiparticles in 
the superconducting state can be obtained from transport measurements, such as thermal 
Hall conductivity^, microwave and terahertz conductivity.^'^ These studies in YBa2Cu307 
(Y123) had revealed a rapid increase of the quasiparticle lifetime on cooling across the super- 
conducting transition, ^'^ and the mean-free-path obtained by the thermal Hall conductivity 
at the lowest temperature attains a value of ~1 //m.^ On the other hand, terahertz con- 
ductivity in Bi2212 shows a scattering rate with rather high values and a slower T-linear 
dependence below Tc, compared to Y123.^ Such differences have been discussed in terms of 
possible impurity effects, which may accompany phase fluctuations, charge and superfluid 
density inhomogeneity, etc. 

Angle-resolved photoemission spectroscopy (ARPES) is a powerful method to investigate 
the momentum-resolved nodal quasiparticle properties, both, as a function of temperature 
and energy [uj). The energy position of the quasiparticle peak in ARPES represents the 
band dispersion, whereas the width indicates the scattering rate P. Recent studies have 
carefully probed the "kink" structure at around 60 meV in the quasiparticle dispersion, 
now commonly observed in hole-doped cuprates.^'^ The kink accompanies a rapid decrease 
of the scattering rate, and is discussed as the renormalization of a bosonic mode, derived 
from phonons or magnons.^'^'^*^ There are important studies^'^^'^^ that discuss the nature 
of the scattering mechanism of the electrons from the uj dependence of the scattering rate, 
e.g. is it a Fermi liquid ? ( with P oc cj^), or a marginal Fermi liquid ? (with P oc |c<j|)^^. 
Due to limitations of momentum- and energy-resolutions, however, the low temperature 
P at o; = (i.e. the Fermi level) has not been obtained accurately enough to discuss 
the nodal quasiparticle properties in terms of unique elementary excitations in the d-wave 
superconducting state. In particular, to date, the discrepancies obtained in its values by 
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transport measurements and ARPES remains to be conclusively clarified. Accordingly, its 
doping dependence at low energy and low temperature is also not yet clarified and remains 
an open question. In this study, we investigate the character of nodal quasiparticles of 
Bi2212 as a function of doping, by VUV laser-excited ARPES measurements. By using a 
low-energy [hv = 6.994 cV) and coherent light source, we can achieve high energy- and 
momentum- resolutions (< 1 meV and 0.0014 A), as well as bulk sensitivity (~ 100 A). 
Prom our measurements, we obtained very sharp quasiparticle peak of 6 meV full-width at 
half-maximum for an overdoped sample. The scattering rate obtained from the quasiparticle 
peak width shows a clear doping dependence, increasing monotonically from the overdoped 
to the underdoped region, reminiscent of the pseudogap phenomenon. It possibly arises from 
the electronic inhomogeneity effect, as is known from scanning tunneling microscopy (STM) 
measurements. 

ARPES measurements were performed using a system constructed with a VG-Scienta 
R4000 electron analyzer and an ultra-violet {hu — 6.994 eV) laser for the incident light. The 
temperature was precisely controlled down to 5 K using a fiow-type He liquid refrigerator. 
The pressure of the chamber was below ~ 5 x 10~^^ Torr throughout all the measurements. 
The energy and angular resolutions were Ej-^g =1.0 meV and ^^j-cs = 0.1°. The Fermi level 
{Ep) of the sample was referred to that of a Au film evaporated on the sample substrate, with 
an accuracy of ± 0.1 meV. High quality single crystals of Bi2212 were grown by traveling- 
solvent fioating-zone method. Tc of the samples from underdoped to overdoped region were 
72 K (UN72K), 80 K (UN80K), 90 K (OP90K), 85 K (OV85K), and 73 K (OV73K). The 
samples were cleaved in situ to obtain clean surfaces. All the data presented in this work 
were obtained at 5 K. 

Figure 1(a) shows the ARPES intensity image along the nodal direction [F - Y: see Fig. 
1(b) for a schematic Brillouin zone] from OP90K. As is well-known,^ the "kink" structure 
refiecting the band renormalization at ~60 meV is clearly observed in the raw data. To 
discuss the detailed features appearing in the band dispersion, we performed a fitting analysis 
of momentum distribution curves (MDC: momentum-profile of the ARPES intensity at a 
fixed binding energy) using the sum of a Lorenzian peak and a constant background, as 
a function of binding energy. The obtained peak positions are overlaid as a red curve in 
Fig. 1(a), and the full- width-half-maximum (FWHM) Ak is plotted as a function of binding 
energy (Eb) in Fig. 1(c), respectively. Ak thus obtained can be approximately described by 
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Ak 2ImE/i>o, where vq and E correspond to the bare Fermi velocity and the self-energy 
of electrons.^ The 60 meV "kink" is easily identified as the knee in the band dispersion. 
In addition, a rapid decrease of Ak at lower binding energies is observed, which has been 
extensively discussed as a renormalization effect due to a bosonic mode. While these "kink" 
properties are very similar to those reported previously by higher-energy ARPES results, the 
absolute value of the linewidth is obviously sharper here ; Ak ~ 0.01 A~^. Similar results 
were recently reported by ARPES experiments using a laser source^^ as well as low-energy 
synchrotron light, demonstrating the effectiveness and impact of low-energy ARPES for 
high resolution measurements, not only in energy but also in momentum. 

In Fig. 2(a-c), the results of the ultrahigh energy- and momentum- resolution measure- 
ment are shown for UN72K, OP90K, and OV73K. A clear decrease of the linewidth on 
increasing the doping level is readily visible from the ARPES images. The energy distribu- 
tion curves (EDCs: energy-spectrum of the ARPES intensity at a fixed momentum) from 
OV73K, shown in Fig. 2(d), indicate the continuous sharpening of the quasiparticle peak 
even at binding energies below £'b = 10 meV. The quasiparticle peak eventually becomes 
a sharp peak with ~ 6 meV FWHM at k = kp [red curve in Fig. 2(d)]. By fitting the 
MDCs by Lorenzians, we obtained clear band dispersions which arc overlaid as red curves 
on the images in Fig. 2(a-c). The FWHM Ak as a function of Eb is plotted in Fig. 2(e) for 
UN72K, OP90K, and OV73K. Although the magnitudes of the FWHM Ak are quite similar 
for OP90K and OV73K (with small differences only at energies below 20 meV) , they are sig- 
nificantly higher for the UN72K sample. In spite of these differences, an overall quasihnear 
LU dependence at low energies is observed for all the samples. Another important observation 
is that Ak at Ep is significantly doping dependent; Ak{0) = 0.014 for UN72K, reduces 
to 0.0062 A-^ for OP90K, and reduces further to 0.0039 A'^ for OV73K. The corresponding 
MDCs at Ep and the fitted Lorenzian curves are shown in Fig. 3 (c). We can estimate the 
mean free path Zmfp of the nodal quasiparticles from Zmfp = Ak{0)~^. The estimated values 
are Z^fp > 260 A (OV73K), 160 A (OP90K), and 70 A (UN72K), respectively. The obtained 
values of Ak and /mfp for all the samples are shown in Table I. It is noted that the results 
for the OP90K sample are in very good agreement with the work of Koralck et al}^ Here 
we would also like to point out that while recent studies have shown the existence of bilayer 
splitting at the nodal point,^^'^^'^^ the MDC widths obtained for OP90K and OV73K in the 
present work are smaller than the reported bilayer splitting of 0.0075 A~^. This indicates 
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that, the present results on OP90K and OV73K, and that of Koralek et al. for optimally 
doped Bi2212, measures only one of the bilayer split bands, presumably the antibonding 
band, at the photon energies of the incident laser used in the studies. The spectral weight 
measured for the bonding and antibonding band is known to strongly depend on the incident 
photon energy. 

Now we compare the EDC at kp among the different doping samples. The EDCs are 
shown in Fig. 3(a). While the EDC from OV73K shows a clear and sharp peak with very 
httle background, the EDC peak from UN72K apparently tends to become broader. To 
get rid of the Fermi-Dirac distribution effect at Ep and estimate the FWHM of the peak 
from EDCs, we symmetrized them at Ep and plotted them in Fig. 3(b). The symmetrized 
EDCs can be well fitted with a single Lorenzian function, as shown by the broken curves 
in Fig. 3(b). The FWHM obtained from the symmetrized EDCs are 12.5±0.6 (UN72K), 
8.4±0.8 (OP90K) and 6.0±0.8 meV (OV73K), respectively, where the errors mainly arise 
from the determination oi kp (see Table I for all samples). The FWHM of an EDC peak 
corresponds to the scattering rate V — h/r ^ vq/vpIvoH, where r and vp correspond to the 
lifetime and the renormalized Fermi velocity of quasiparticles. The relation between EDC 
and MDC peak FWHM can be expressed as F a; vpAk, where vp can be estimated from 
the gradient of the band dispersion in the vicinity of Ep. By linearly approximating the 
dispersion at lOmeV < -Br < meV obtained from the MDC analysis (red curves in Fig. 
2(a-c)), we can get vp within errors of about ±10% and thus F from MDCs in Fig. 3(c) as 
15±1.4 (UN72K), 8.1±0.6 (OP90K), and 5.1±0.4 meV (OV73K), which are very similar to 
the values obtained from EDCs (see Table I). Thus, we succeeded in obtaining the coherent 
nodal quasiparticle component from ARPES and estimating its scattering rate F. Here, it 
is important to compare our results with that from previous transport measurements. The 
quasiparticle scattering rate obtained from terahertz conductivity using a nearly optimally- 
doped sample {T^ — 85 K)^ is 1/tqp ~ 3 THz ~ 12 meV at the lowest temperature of 
measurement, 12 K. By assuming a simple linear extrapolation, the scattering rate should 
be I/tqp ~ 2 THz ~ 8 meV at around 5 K. This value is equivalent with our result from 
OP90K (8.1 - 8.4 meV), indicating full consistency among two completely different probes. 
While not yet reported, a similar doping dependence of 1/tqp can also be expected in 
terahertz conductivity measurements on Bi2212. 

The doping dependence of the FWHM F thus obtained is plotted in Fig. 4, with that 
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of Tc- r shows a monotonic increase from overdoping to underdoping. Such tendency it- 
self is very similar to that of the well-known pseudogap in (7r,0) region, being intensively 
discussed to date as the manifestation of antiferromagnetic spin fluctuation, charge order, 
phase fluctuation, and so on. Very recently, systematic investigations of the electronic struc- 
tures by Raman scattering^^ and ARPES^^'^^ measurements have revealed that there is a 
strong dichotomy among the near-nodal (Fermi arc) and the antinodal (pseudogap) regions. 
The Fermi arc region shows the opening of a well-defined d-wave superconducting gap with 
coherent quasiparticles regardless of the doping level, while the pseudogap shows a mono- 
tonic increase on underdoping with the enhancement of the incoherent spectral shape. STM 
measurements also indicate the spatially inhomogeneous (pseudo)gap spectra arising in the 
underdoped samples, while keeping rather homogeneous coherent quasiparticlc spectra for 
all doping levels. ^^'^^'^^ Our observation of the well-defined coherent peaks is basically in 
agreement with this picture. However, it also shows that the lifetime of the nodal quasipar- 
ticlc at the center of the Fermi arc shows a significant doping dependence, when investigated 
with high energy and momentum resolution at T = 5 K. 

The lifetime of the nodal quasiparticle can be affected by a number of scattering mecha- 
nisms. In case of a metal with nearly flat density of states (DOS) near Ep, electron-electron 
scattering can provide the low-energy scattering rate (a; 0) as ImE oc max(ct;^,T^) in a 
Fermi liquid, while the ImE oc max{\u!\,T) in a marginal Fermi liquid picture.^^ Similarly, 
the electron-phonon scattering gives ImE oc max(a;^,T^). In a d-wave superconducting 
state, they are completely modified refiecting the |a;|-linear DOS near Ep."^^ The dissipation 
by electron-electron (or spin fluctuation) scattering is given by ImE oc max(ct;^, T^), and 
electron-phonon by ImS oc max(a;^,T^). In the low energy and temperature limit (appro- 
priate for the present measurements performed at T < 0.07Tc with Ej-cs ~ 1 meV), however, 
all the above scattering effects should become negligible. In fact, the cu dependence of the 
ImE oc Ak shows quasilinear behavior, indicating that the above dissipation mechanisms 
are not dominant in this T,u! region. One of the remaining terms is the elastic impurity 
scattering mechanism, which can affect a c?-wave superconducting state in a non-trivial way, 
unlike in a normal metal. ^'^^ An in-plane impurity, for example, is known to strongly scatter 
the quasiparticles and create a localized state around the node appearing as an enhancement 
of ImE (a;) near Ep.^^ In case of the out-of-plane impurity which likely provides a weak for- 
ward scattering potential, on the other hand, ImE of the quasiparticles can show a |a;|-linear 
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dependence at a; — > with the value of ImE(a; = 0) depending on the scattering length of 
impurity potential and its concentration.^^ Ak showing a cu dependence of convex upward 
structure at around = 10 meV in OP90K may be due to such an elastic scattering, as 
expected from calculations based on a BCS rf-wave model. Recently, the interplay among 
the apical oxygen dopant, the gap inhomogeneity and quasiparticle interference has been ob- 
served in STM measurements. The weakening of the screening effect due to the reduction 
of carrier density may make the system more susceptible to impurity scattering, typically 
emerging as granular (typical size of 3 nm) superconductivity in underdoped samples^^'^^. 
It is interesting that the mean free path of the nodal quasiparticle in UN72K is l^fp — 70 A, 
a little bit longer than the reported size of the granular superconducting domains. On the 
other hand, we cannot yet completely rule out the possibility of the quantum fluctuation 
enhancement in underdoped samples, such as phase fluctuation arising from the reduction 
of superfluid density, charge order fluctuation, etc, if they have very low energy scales (< 
1 meV). Future work including detailed temperature-dependent measurements are required 
to elucidate on these issues. 

In conclusion, we have performed an ultrahigh resolution angle-resolved photoemission 
spectroscopy measurement to elucidate the properties of nodal quasiparticles in a rf-wave 
superconductor. The quasilinear cu dependence of the scattering rate F obtained from MDCs 
suggest the role of an elastic scattering mechanism by out-of-plane impurities, such as oxygen 
dopants. The value of F at a; = for optimally and overdoped samples are 8 meV and 
6 meV, respectively, in a good correspondence with terahertz conductivity results. The 
doping dependence of F resembles that of the pscudogap, monotonically increasing up to 12 
meV on approaching the underdoped region. The estimated mean free paths accordingly 
decreases from 260 A in an overdoped to 70 A in an underdoped sample. This suggests 
role of electronic inhomogeneity on the nodal quasiparticle scattering at low temperature in 
underdoped samples, as is known from scanning tunneling microscopy (STM) measurements. 

We thank Y. Yanase for valuable discussion. 
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sample 


P 


Afc [A-i] 


^mfp [A] 


VF [eVA] 


vp^k [meV] 


r [meV] 


UN72K 


0.11 


0.0136 


74 


1.1 (±0.1) 


15.0 (±1.4) 


12.5 (±0.6) 


UN80K 


0.12 


0.0136 


74 


1.1 (±0.1) 


15.0 (±1.4) 


11.7 (±0.6) 


OP90K 


0.16 


0.0062 


160 


1.3 (±0.1) 


8.1 (±0.6) 


8.4 (±0.8) 


OV85K 


0.19 


0.0075 


130 


1.0 (±0.1) 


8.0 (±0.8) 


7.2 (±0.6) 


OV73K 


0.21 


0.0039 


260 


1.3 (±0.1) 


5.1 (±0.4) 


6.0 (±0.8) 



TABLE I: Doping dependence of the nodal quasiparticle property at a; = 0. p is the hole concen- 
tration estimated from % = r™^(l - 82.6(0.16 -pf), Ak, /^fp, vp, and vpAk are FWHM, mean 
free path, Fermi velocity, and scattering rate obtained from MDC analysis, T is the scattering rate 
obtained from EDC analysis. See the text for details. 

FIG. 1: (color), (a) ARPES intensity image plot from Bi2212 OP90K. The red curve shows the 
band dispersion as obtained from the MDC peak positions, (b) A schematic of the Brillouin zone 
of Bi2212. The red line shows the momentum-region of the measurement [along (0,0)-(7r,7r)]. (c) 
£;B-dependence of the MDC FWHM obtained from OP90K. 

FIG. 2: (color), (a)-(c) Ultrahigh resolution ARPES intensity image plot from Bi2212 UN72K, 
OP90K, and OV73K, respectively. The red curves show the band dispersions obtained from MDC 
peak positions, (d) EDG of taken from ARPES on OV73K (c). (e) £;B-dependence of MDC FWHM 
obtained from UN72K (a), OP90K (b), and OV73K (c), respectively. 

FIG. 3: (color onhne). EDC at k = kp (a) and that symmetrized at Ep (b) from ARPES on 
Bi2212 UN72K, OP90K, and OV73K [Fig. 2(a)-(c)]. (c) shows the MDC at E = Ep obtained 
from the same ARPES spectra. The curves in (b) and (c) are the fitting result by using a single 
Lorenzian function. 
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FIG. 4: (color online). Doping dependence of the scattering rate T estimated from the FWHM of 
EDO and MDC in Bi2212. The right axis shows the of the measured samples with the relation 
to hole concentration p, = T™'^(1 - 82.6(0.16 - p)^). 



11 




Figure 1 (K. Ishizaka et al) 
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